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Three-dimensonal pentamode metamaterials are artificial solids that effec- 
tively behave much like liquids, which have vanishing shear modulus. Pen- 
tamodes have recently become experimental reality. Here, we calculate their 
phonon band structures for various parameters. Consistent with static con- 
tinuum mechanics, we find that compression and shear waves exhibit phase 
velocities that can realistically be different by more than one order of mag- 
nitude. Interestingly, we also find frequency intervals with more than two 
octaves bandwidth in which pure single-mode behavior is obtained. Herein, 
exclusively compression waves exist due to a complete three-dimensional band 
gap for shear waves and, hence, no coupling to shear modes is possible. 
Such single-mode behavior might, e.g., be interesting for transformation- 
elastodynamics architectures. 

Artificial materials called metamaterials have recently started to emerge in acoustics 
and elasto dynamics, for example see [IJ[2J[3lEl[5j[6l|71[8j|9]. This field is driven by 
the quest for more control on elastic/acoustic waves, inspired by corresponding progress 
on three-dimensional metamaterials in electrodynamics/optics (for recent reviews see 
|10| 111)). Three-dimensional (3D) pentamode metamaterials can be seen as the elasto- 
dynamic counterpart of 3D magnetodielectric metamaterials in optics in the following 
sense: In 3D transformation optics |12[ I13j. magnetodielectrics enable anisotropies with- 
out birefringence and, hence, an optical wave impinging onto some structure with a 
certain polarization also emerges with that polarization and does not generate partial 
waves with other polarizations. In mechanics or, more precisely, in elastodynamics, nor- 
mal spatially inhomogeneous elastic solids would couple compression and shear waves 




Figure 1: (a) Illustration of the structure approximating the pentamode metamaterial 
ideal considered in this Letter. Both the extended face-centered-cubic (fee) unit 
cell (red) and the primitive cell (black) are shown, (b) Corresponding body- 
centered-cubic (bec) Brillouin zone. The characteristic points are indicated, 
aiming at easing the interpretation of the band structures in Fig. 2 and Fig. 4. 

[U HH ESI HH [T71 HE]. Thus, an incident compression wave would lead to emerging shear 
waves as well, making the realization of transformation-elastodynamic structures in three 
dimensions impossible. Pentamodes are a mathematical ideal |14} H9] for which the shear 
modulus is strictly zero, thus no shear modes exist at finite frequencies. Recently, we 
have fabricated 3D microstructures [7] that approach the pentamode ideal in that the 
effective bulk modulus computed within static continuum mechanics is more than three 
orders of magnitude larger than the effective shear modulus. 

However, it is not obvious how these static results can be translated to the dynamic 
case, i.e. , to phonon waves. For example, one would like to know the phase velocities of the 
different modes. The moduli are known from static continuum mechanics. For obtaining 
knowledge on phase velocities, one additionally needs the effective mass densities. It 
is known that the dynamic effective mass density can be distinct from the static mass 
density as one would determine by using a mass balance (see, e.g., |20[ [2TJ I22j). 

Furthermore, it has not been clear up to which maximum operation frequency penta- 
modes would keep their desired effective material properties. Band structures have been 



2 



published for honeycomb-lattice structures |23| that may be seen as the two-dimensional 
counterpart of pentamodes (i.e., as bimodes), but not for three-dimensional pentamodes 
so far. Thus, in this Letter, we compute and discuss pentamode phonon band structures 
systematically as a function of the underlying structure parameters. To connect to our 
previous static work [7] and to allow for a direct comparison, we choose parameters for 
the constituent material corresponding to our recent work. 

Fig. [ija) recalls the structure approximating the pentamode ideal suggested by Milton 
and Cherkaev |19| . Both the extended face-centered-cubic (fee) cell with lattice constant 
a and the primitive cell (black) are depicted. Cones touch each other at their thin ends 
with diameter d. These connection points form a diamond lattice. Pairs of cones are 
connected at their thick ends with diameter D. The length of the double-cones is then 
given by h = \/3a/4. Panel (b) of Fig. [I] illustrates the corresponding body-centered- 
cubic (bec) Brillouin zone (i.e., that of diamond). The usual characteristic points are 
indicated, aiming at easing the interpretation of the corresponding band structures shown 
in Fig. [2] and Fig. |] 

For the numerical calculation of the pentamode band structure, we choose a polymer as 
constituent material with Young's modulus 3 GPa, Poisson's ratio 0.4, and mass density 
p = 1190 kg/m 3 (as previously [7J). We solve the elastodynamic equations |14j in vacuum. 
The results should also be applicable to pentamodes in air as long as the air speed of 
sound does not match that of the pentamode branches, in which case strong coupling 
could occur. The coupled elastodynamics-acoustics problem has only been addressed 
in very rare cases |24j . The scalability of the elastodynamic equations allows to easily 
translate our results to other choices of the Young's modulus, provided the Poisson's ratio 
as well as the structure parameters are fixed. For example, if one increases the Young's 
modulus by factor 100, the frequency increases by the square root, i.e., by factor 10. We 
use a commercial software package (COMSOL Multiphysics, MUMPS solver) with Bloch 
boundary conditions imposed for the primitive real-space cell shown in Fig. [lja). We 
have carefully checked that all results depicted are converged. Typically, convergence is 
achieved by using around 4 x 10 4 tetrahedra in one primitive real-space cell. Also, we 
choose a = 37.3 \xm throughout as previously [7]. The right-hand side vertical frequency 
scale in Fig. [2] is in S.I. units. On the left-hand side vertical scale, we also show the a/A 
ratio as a normalized and more universal frequency scale. Here, A is the wavelength of 
sound in air, for which we have taken a standard sound velocity of 343 m/s. For example, 
a ratio of a/A = 0.1 means that the air wavelength is ten times larger than the lattice 
constant a (whatever the value of a might be). The usual tour through the Brillouin 
zone on the horizontal axis is parameter-independent anyway. 

In Fig. [2j we fix D = 3 (J.m and vary d from (a) 3 yon. to (b) 0.55 um to (c) 0.2 \xm 
(compare insets in Fig. [3]). Note the different vertical scales. The branches emerging from 
the r point (i.e., wave vector = (0, 0, 0), see Fig. [l|h)) are the dispersions of the "acous- 
tic" shear and compression waves respectively. Their slopes are the corresponding phase 
velocities (see red straight lines). Obviously, the ratio of the slopes increases with de- 
creasing diameter d. This trend was to be expected from the static continuum-mechanics 
calculations in Ref. [TJ. A large ratio of the slopes is desirable for transformation elas- 
todynamics, because it suppresses undesired coupling between shear and compression 
waves in spatially inhomogeneous architectures. We will come back to the dependence of 
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Figure 2: Calculated band structures for the pentamodes shown in Fig. |1] The fixed 
parameters are a = 37.3 [im (hence h = 16.15 |xm) and D = 3 (Jin. The 
diameter of the connection region d (see Fig. [TJ is varied from (a) 3 \im to 
(b) 0.55 (J.m to (c) 0.2 [im. The red straight lines are fits to the dispersion 
branches of interest in this Letter. They correspond to phase velocities of (a) 
cg = 198 m/s and cb = 596 m/s, (b) cq = 25m/s and eg = 363 m/s, and (c) 
cq = 6 m/s and cb = 240 m/s. In the regions highlighted in gray, only a single 
phonon mode occurs due to a complete three-dimensional band gap for shear 
waves. By the scalability of the problem one can, for example, simultaneously 
replace \xm by mm and MHz by kHz. 
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Figure 3: The ratio cb/cq (red) of the phase velocities of compression (black) and shear 
waves (compare red straight lines in Fig. 2) for different parameter combina- 
tions of d and D (see Fig. [TJ as derived from band structure calculations like 
shown in Fig. [2] for the TX direction (i.e., propagation of the waves along any 
of the six face diagonals of the cube in FigQa)). The volume filling fractions 
/ are given in green. 

this ratio below. Furthermore, the gray regions in Fig. [2] highlight frequency intervals in 
which only a single phonon mode exists in the pentamode structure. The relative width 
of the lower frequency interval increases with decreasing diameter d from (a) to (c) in 
Fig. [2] In other words, it increases as one approaches the pentamode ideal. Our previous 
experiments j7] have realized about d = 0.55 ]xm. In the corresponding Fig. |2^b), the 
upper frequency bound of this interval is more than two octaves larger than its lower 
bound. In these single-mode regimes, only compression waves exist and shear waves are 
suppressed by a complete three-dimensional band gap. We will speculate about possible 
applications of this finding at the end of this Letter. At yet higher frequencies, we find 
the usual "spaghetti" region. Many rather flat bands describe more spatially localized 
modes that are not of particular interest to us in the context of this paper. 

The dependence on the large diameter D is much less pronounced. To give an overview, 
the matrix structure in Fig. [3] shows the ratio of the phase velocities of compression and 
shear waves for various combinations of D and d. In each box of this matrix, the structure 
is illustrated too. The phase- velocity ratio is nearly independent on D and scales roughly 
inversely with d. This qualitative scaling is consistent with our previous static continuum- 
mechanics calculations [7j. There, we found that the ratio of the bulk modulus B and 
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the shear modulus G scales approximately like 




In continuum mechanics, the phase velocity of the compression wave, cr, is given by 
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with the effective (dynamic) mass density In the last step, we have assumed B/G 3> 
1. For the shear waves, one has 

CG = 

Naively, one would assume that both mass densities Pb = PG = Po can be taken as the 
static mass density po, which is what one would measure with a mass balance in the 
laboratory. Mathematically, po = fp is simply given by the volume filling fraction / 
times the mass density p of the constituent material. The more sophisticated Berryman 
formula [20j is often a more accurate approach but is also known to fail quantitatively and 
qualitatively in some cases |21j . In any case, when taking the ratio cb/cg, any common 
mass density would drop out and one gets a scaling according to cb/cg oc 1/(cLD 1//4 ) 
oc 1/d. Indeed, we obtain a good least-squares fit to the numerically computed values 
in Fig. 3 with the formula 

cr h 
-RfO.5-: . 
cg d 

For example, for d = 0.55 |im, D = 3 |0.m, and a = 37.3 |J.m, hence h = 16.15 |xm, we have 
cb/cg ~ 16 in Fig. [3] For the same parameters, the above heuristic continuum-mechanics 
formula yields B/G ~ 1260 [7J, which, assuming identical mass densities ps = PGi leads 
to cb/cg ~ v 1260 ~ 35. Considering the simplicity of this reasoning, the result is not 
too far off. One should also be aware that there is no unique phase velocity of the shear 
waves as can be seen from the splitting of the corresponding dispersion branches in the 
rK direction in Fig. [2] We will come back to the aspect of isotropy below. 

In static continuum mechanics [7], we found very little dependence of the results on 
the Poisson's ratio of the constituent material (the polymer). We have also investigated 
this aspect for the band structures. For example, when drastically changing the Pois- 
son's ratio from 0.4 (as above) to 0.1 and keeping all other parameters fixed, the relevant 
frequencies in the band structure shift by less than 3 % (not depicted). However, chang- 
ing the Young's modulus of the constituent material will, of course, stretch/compress 
the frequency scale by a common factor as discussed above but it will not change the 
qualitative behavior of the band structure at all. Likewise, a change in lattice constant 
a leads to a trivial inverse scaling of the frequency axis. 

Let us come back to the above mentioned single-mode regime. In Fig. Qa) for d = 
0.55 (J.m, its upper limit lies at a normalized frequency of about a/A = 0.4 (the lower limit 
lies below a/A = 0.1). This means that the corresponding wavelength in air is 2.5 times 
larger than the pentamode fee lattice constant a at the upper limit. Below this upper 
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Figure 4: Pentamode dispersion relations for four selected characteristic directions in the 
first Brillouin zone (compare Fig. [l|h))- The parameters and the gray single- 
mode region are as in Fig. \2\ 



frequency, the dispersion can still be approximated reasonably well by a straight line. 
One might expect that the effective medium approximation is still reasonably well fulfilled 
here too. Also, we would still expect a more or less isotropic behavior of the compression 
waves - as for an ideal pentamode |19| . To investigate this aspect, we depict in Fig. [4] the 
numerically calculated dispersion relations for various directions in momentum space for 
d = 0.55 (J.m (i.e., all parameters as in Fig. (2^b)). Obviously, the steeper compression- 
wave branches in Fig. [4] exhibit only very little dependence on direction, i.e., the behavior 
is very nearly isotropic in 3D. In contrast, differences in the slopes as large as a factor of 
two occur for the flatter shear-wave branches. Indeed, it was already apparent from the 
lifting of degeneracy in the dispersions in the TK direction in Fig. [2] that the pentamode 
shear-wave dispersion is not isotropic at all. 

Finally, we speculate that the single- mode region (Fig. [2]) with nearly isotropic behavior 
for the compression waves and with a complete three-dimensional band gap for shear 
waves might be interesting for future three-dimensional transformation-elastodynamic 
architectures. After all, one motivation for constructing pentamode metamaterials lies 
in suppressing the shear waves in 3D. One should be aware though that the different 
local single-mode regions of a generally complex, spatially inhomogeneous and possibly 
also anisotropic pentamode-based elastodynamic architecture will likely not fully overlap. 
However, one could choose the operation frequency to be within the single-mode region 
in the outer parts of the overall structure. In this case, at least no shear waves could 
emerge, e.g., from a cloak. 

In conclusion, we have calculated phonon band structures for 3D metamaterials ap- 
proximating the pentamode ideal for the first time. Our findings agree qualitatively with 
our previous static continuum-mechanics calculations in that we find ratios of the phase 
velocities of compression and shear waves as large as 16 for experimentally accessible 
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structure parameters. Moreover, we also find an interesting single- mode frequency in- 
terval with a width larger than two octaves. For reasonable material parameters, its 
upper frequency bound corresponds to an air wavelength that is only 2.5 times larger 
than the pentamode fee lattice constant. This means, for example, that order a = 1 cm 
pentamode lattice constants would suffice for up to 10 kHz operation frequencies. Within 
this single-mode regime, the compression-wave dispersion is still nearly isotropic in three 
dimensions. 

We thank Andreas Frolich, Jonathan Miiller, and Michael Thiel for help with imple- 
menting the periodic boundary conditions. We acknowledge support by the DFG-Center 
for Functional Nanostructures (CFN) via subproject A1.5. 
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